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Abstract 

Mesoscopic particles immersed in a critical fluid experience long-range 
Casimir forces due to critical fluctuations. Using field theoretical methods, 
we investigate the Casimir interaction between two spherical particles and 
between a single particle and a planar boundary of the fluid. We exploit the 
conformal symmetry at the critical point to map both cases onto a highly 
symmetric geometry where the fluid is bounded by two concentric spheres 
with radii R- and In this geometry the singular part of the free energy 
only depends upon the ratio and the stress tensor, which we use 

to calculate has a particularly simple form. Different boundary conditions 
(surface universality classes) are considered, which either break or preserve 
the order-parameter symmetry. We also consider profiles of thermodynamic 
densities in the presence of two spheres. Explicit results are presented for 
an ordinary critical point to leading order in e = 4 — d and, in the case of 
preserved symmetry, for the Gaussian model in arbitrary spatial dimension d. 
Fundamental short-distance properties, such as profile behavior near a surface 
or the behavior if a sphere has a 'small' radius, are discussed and verified. The 
relevance for colloidal solutions is pointed out. 

PACS numbers: 05.70.Jk, 68.35.Rh, 75.40.Cx, 75.30.Pd, 82.70.Dd 
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I. INTRODUCTION 



Consider a fluid system close to a critical point, for example a fluid near the liquid- 
gas critical point, a binary liquid near the consolute point, or liquid ^He near the A-point. 
Due to the long correlation length a local perturbation, as provided by the wall of the 
container or by the surfaces of immersed colloid particles, is not screened within a few 
atomic layers but influences the system over appreciable distances. This gives rise to a 
long-range interaction between immersed particles or a particle and a wall, as pointed out 
by Fisher and de Gennes |l|J^. The forces may be termed critical Casimir forces [§] on 
account of the close formal relationship with the usual Casimir forces |^ due to quantum 
fluctuations of the electromagnetic field. An interesting feature of the critical Casimir effect 
is that it can be turned on and off by changing the magnitude of the correlation length, e.g. 
by changing the temperature. 

The simplest and most-considered geometry for the critical (as well as the electromag- 
netic) Casimir effect is the parallel-plate geometry. Results have been obtained in spatial 
dimensions c? = 2 by using conformal invariance and in c? = 4 — e |ll0|-|l3[]. In this paper 



we consider spherical boundary surfaces, i.e. spheres immersed in a critical fluid. Many of 
the general model-independent properties for this geometry are also discussed in a recent 
letter by Burkhardt and one of the present authors . 



The spherical shape is of immediate relevance for colloidal particles |jT5|,[T6|, and besides 
the planar boundary it is the simplest boundary shape for theoretical analyses. The reason 
is that systems right at the critical point display an invariance not only with respect to 
homogenous dilatations but also with respect to conformal or angle-preserving coordinate 
transformations [llTi. The class of special conformal transformations ||18| which exists in 



arbitrary spatial dimension d maps spheres onto spheres. This has been used to study the 
critical properties of an infinite system with a single spherical hole (particle) and of a finite 



system with a spherical boundary by conformally mapping results for the half-space |ll9| , |20 

An infinite system with two arbitrary non-overlapping spherical holes (particles) may be 
obtained by conformally mapping a highly symmetric finite system which is bounded by 
two concentric spheres |2T[]. This is shown in Fig. 1 and follows from the special conformal 
transformation 

r' r + K R 

— = ^ — (1.1) 

r'2 |r + R|2 2i?2 ^ ^ 

which maps the sphere of radius R centered at the origin of r-space onto a plane in r'-space 
that is perpendicular to R and passes through the origin. Two spheres with fixed radii R-, 
centered about the origin of r-space (and depicted in Fig. 1 a) map onto two spherical 
particles in r'-space, as shown in Fig. 1 b. The points of intersection I, II, III, IV with the 
r^-axis have the coordinates 

f ' ' ' ' \ - oD ( R+ + R R+- R R- R- R + R^ \ 
{r^i,r^,,,r^,,,,r^,^)-lKy + ^' ^ + ^_ ' ^ _ ^ J l^-^i 

with R_ < R< R+. 

Two cases of particular interest will be studied in detail below: 
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(i) For R = y/R_R+, the two spheres in r '-space have equal size (SES), since r^j = —r'j_iY 
and riii = -rini- 

(ii) For R — > the left sphere reduces to a plane through the origin, which is denoted 
by a broken line in Fig. lb. This corresponds to one sphere with radius -Ri near a 
planar wall (SPW) of the container. 

One may also consider a single spherical particle immersed in a critical fluid inside a 
spherical container. This geometry is obtained from the concentric geometry via ( [I.ID when 
R is chosen outside the interval For R > R^ the order of the intersection points 

is II, III, IV, I, as in the concentric geometry of Fig. 1 a, and sphere a remains inside sphere 
b. For < -R < R- the order is IV, I, II, III, and sphere a is located outside sphere b. 

An important parameter of the two-sphere geometry is the cross ratio of intersection 
points, which remains invariant under conformal transformations [0,0 . Instead of working 
directly with the cross ratio, it is convenient to introduce the related invariant quantity 

^=i^z /f^ t w" (1-3) 



yi.iv '^±iiiJ V ±11 ^±1) 



which can be written as 



-= ^'\^'j, ^' (1-4) 

in either of the above geometries. Here Ri and R2 are the radii of spheres a and 6, re- 
spectively, and S12 is the distance between their centers. On comparing with the concentric 
geometry of Fig. 1 a, where S12 = and where Ri and R2 have been denoted by i?_ and 
respectively, one finds from (ITT 



with 



'^ = ^(p + pO' 1</«<oo, (1.5) 



P=^, 0<p<l. (1.6) 

One can also obtain ( |1.5| ) by inserting ( |1.2| ) into ( |1.3| ), which demonstrates the conformal 
invariance of k since the parameter R of (|1 . 1| ) drops out. It follows from ( |1.3|) or 



that K = 1 corresponds to two particles that touch, while k — > 00 describes the situation 
where one or both of the radii R2 are much smaller than the remaining lengths in the 
two-sphere configuration. 

Now let us turn our attention to the free energy. For two non-overlapping spherical 
particles in an infinite critical medium we consider 

(5F(S12, R2 ) = Fa,tiSl2, Rl, R2 ) - UaiRl) - f/fe(i?2) , ^12 > Ri + R2 ■ (1.7) 



Here F is the singular part of the free energy for immersing two particles into an 
infinitely extended fluid, while U is that for immersing a single particle. 6F determines 
the Casimir force, —d 6F/dsi2, and it vanishes as the separation Si2 tends to infinity with 
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R2 fixed. For the single particle (a, Ri) in the fluid with spherical confinement (6, -R2) 
one may consider a similar expression (5F(si2, -Ri, R2) with S12 + -Ri < -R2, where Ua{Ri) 
and Ub{R2), the free energy of the spherical container without a particle, are subtracted 
from the total free energy. This expression vanishes as R2 tends to infinity with Ri and S12 
fixed. Right at the bulk critical point the quantities 6F are invariant under special conformal 
transformations [0,^ and can only depend on the conformal invariant n in ( |1.4|) : 



6F{si2, Ri, R2) = kBT,J^{K) (1.8) 

It follows from the above arguments that 6F for two particles in the infinite critical liquid 
and 6F for the single particle in the spherically confined liquid both obey (|1.8|) with the 
same function We note two special cases. First, for a single sphere with radius Ri and 

with its closest point at a distance D from a planar wall (SPW), the Casimir interaction is 

6Fsp^{D, = fc^T, ^ (1 + D/Ri) . (1.9) 

This follows from (|1.8|) by substituting S12 = Ri + R2 + D and with i?2 — ^ 00. Second, for 
the concentric geometry (CON) of Fig. 1 a, the result reads 

SFco^iR_,R+) = kBT,J^{K{R./R+)) , (1.10) 



with k(p) given by (|1.5|) . Due to the high symmetry, the concentric geometry is especially 
convenient for the explicit evaluation of the universal scaling function JF(k). 

The scaling function JF depends on both the bulk universality class and the surface 
universality classes a, h pB| , [2B[ ] of the critical system. For a liquid-gas critical point or a 
critical consolute point of a binary mixture, a and h correspond to the Ising model with sur- 
face magnetic fields [P5|-PB| either parallel (t, t) or antiparallel (T, i) depending on whether 



the two surfaces favor the same or different signs of the order parameter. Liquid "^He at the 
A-point belongs to the bulk universality class of the XY model, and surfaces are expected 
to suppress the superfiuid order |]29|, corresponding to the surface universality class O of 



the 'ordinary transition' ||25|j26[l . One may also consider tricritical bulk universality classes 
both for mixtures of ordinary liquids and of ^He-^He. It is also interesting to consider the 



multicritical surface universality class SB of the 'special transition' P5| , P6 |. 

Exact results for JF have been obtained by Burkhardt and Eisenriegler [|I^ for the critical 
Ising model with pairs of surface universality classes (a, h) = (O, O), (O, t), (T? T); and (|, |) 
in spatial dimension d = 2. In the present paper we investigate JF for d near the upper 
critical dimension dc = 4. For the combinations (a, 6) = (T,?), (T,!), ihO), and {^,SB), 
where the order parameter symmetry is broken, the leading behavior in the e-expansion 
follows from mean-field theory. This will be considered in Sec. Here we make close 



contact to the work of Indekeu et al. |]rT| for the parallel-plate geometry and of Gnutzmann 
and Ritschel for the order-parameter profile in the concentric geometry for (t, T)- For 
the symmetry-preserving combinations (O, O), {O, SB), and {SB, SB), the leading behavior 
is that of the Gaussian model. Since the Gaussian model is conformally invariant for all 
d > 2, it is instructive to consider the above pairs of boundaries in general d. This is carried 
out in Sec. |IV|. The parallel-plate result (Krech and Dietrich |]T2[) emerges as a special case. 
The explicit evaluation of the different universal scaling functions JF in Sees. |T| and 



IV| are preceded in Sec. |T| by a discussion of general relationships involving JF which hold 
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for all spatial dimensions and universality classes of interest. We discuss the limits k — > 1 



and K — CX3, which have also been addressed in Ref. |]14[ and correspond to p — 1 and 
p ^ 0, respectively. For p — > 1, the scaling functions are determined by the amplitudes 
of the Casimir effect in parallel-plate geometry [p|-[l3[|. For p ^ 0, we use the small-sphere 
expansion ||T^ (a kind of short-distance expansion) to relate JF to the amplitudes of the bulk 
correlation function and of the profile in the half-space of either the order parameter or the 
energy density, depending on the pair a, b of universality classes. We also discuss how JF can 
be conveniently calculated from the thermal average of the stress tensor Furthermore, 
we discuss the profiles of thermodynamic densities in the presence of two spheres |^T|. For 
the particular case of a sphere and a planar wall (SPW), we use a short- distance expansion 
ll30| , |3llJT3| to relate the behavior of the profile near the wall to the average of the stress tensor 
at the wall, i.e. to JF. All the general relationships are consistent with the explicit results 
in Sees. |ITl| and |V. 

Sec. |V| contains our conclusions. 



II. GENERAL RELATIONSHIPS 

Here we discuss several important relationships which involve the scaling function J-' of 
the free energy. They are of a general character and apply to most of the bulk and surface 
universality classes and spatial dimensions d which are of interest. 



A. Free energy and the stress tensor 

While conformal coordinate transformations of the two-sphere geometry preserve k (and 
thus leave 6F invariant), certain non-conformal transformations preserve the spherical shape 
of the boundaries but change k. In equation ( |A^) we consider a transformation which 
increases the distance su of the two non-overlapping spheres of Fig. 1 b keeping Ri and R2 
fixed. This is directly related to the Casimir force. In the concentric geometry of Fig. 1 a 
one may apply the infinitesimal transformation r ^ f defined by 

f — r = ar6(r — ro) , a <^ 1 (2.1) 

which increases the radial component of r by a fraction a if r > tq. Choosing ro between 
and one sees that increases while i?_ remains unchanged. Consequently, p and 
K are changed. The effect of these transformations on 6F can be written as an integral over 



the stress-energy tensor Tm [0. In the concentric case 

P 6FcoM = ksTc Sd r^ (T„„(ro))coN , (2.2) 
dp 

as we show in more detail in Appendix A. Here 5*^ = 2%'^^'^ /T{d/2) is the surface area of 
the (i-dimensional unit sphere, and T„„(ro) is the radial component (normal to the spherical 
surfaces) of the stress tensor at a distance tq from the center. The thermal average at the 
right hand side of Eq. ( |2.2| ) is taken at the critical fixed point, where the trace of the stress 
tensor vanishes |T^. Since ro is arbitrary in the interval (i?_,i?+). 
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(T„„(ro))coN = (2ro)-'^y(p) , 



(2.3) 



where we have included a factor 2"'^ for convenience below. In evaluating J-'{k) and the 
Casimir forces in Sec. |T| and we have found it convenient to first calculate the scaling 
function Y defined in (p.3|) and then integrate the result according to (p.2|). 



B. Free energy for k ^ 1 

The limit p 1 corresponds to the parallel-plate geometry. The concentric geometry 
becomes equivalent to the parallel-plate geometry if we let i?+ —>■ oo while keeping the width 
L = i?+ — fixed. In this limit 1 — p = L/R^ tends to zero. With the usual notation 
0JTl|Jl^ Aa^b/L'^~^ for the Casimir free energy per fc^T and surface area of one of the two 



parallel plates a, b, 

SFaoM kBT, Aa,t (1 - p)'-' Sd (2.4a) 

and 

F(p)^A,,,(f/-l)2'^(l-p)-'^ (2.4b) 

as p ^ 1. Eq. (|2.4a|) also determines the Casimir interaction of non-overlapping spheres that 
nearly touch. This follows from the form SdAa,b [2(/t — l)]"'-'^"^-'/^ of JF for 0<k — l^las 
implied by ( |1.1CI| ) and ( |1.5| ). In this geometry the limiting form reads [|T^] 



5F{si2, Ru R2) ^ kBT, Sd Aa,b \2iR^' + R^') d] . (2.5) 



Here (|L^) and (|OD have been used, D = su — Ri — R2 is the distance between the 
closest points of the two spherical surfaces (i.e. r_J_jjj — r^jj in Fig. lb), and we assume 
D Ri, R2. In the limit i?^^ 0, the right hand side of (p.5|) determines the Casimir 
interaction (5Fspw(-D, -Ri) of a single sphere close to a planar wall {D <^ Ri). 



K ^ 00 



C. Small-sphere expansion and free energy for 
In the other limit p — 0, where the radius of the inner sphere is much smaller than 



that of the outer boundary, one may use a short-distance expansion introduced in Ref. |]T4 
In this expansion the Boltzmann weight exp(— 7^5) which generates in the critical system 
a spherical hole S with center at and radius Rs is written in terms of local fluctuating 
quantities (operators) S(r5) and non-fiuctuating amplitudes X in the form 



1 + a (2.6a) 



(e-^5)buik 
with 

a = J2 MRs) H(r5) . (2.6b) 
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The relation ( |2.6| ) holds in thermal averages or correlation functions containing other fluctu- 



ating quantities (operators) or boundaries [32| with distances from that are much larger 
than Rs- We mention a few examples where ( p.6| ) is useful: Consider the thermal average 
in bulk of a primary operator \E'(r) such as the order parameter or the energy density 



in an inflnite critical system (bulk) with a spherical hole S. By means of ( p.6a|) one finds 



(^(r)) 



5 in bulk 



(e-^« ^(r))buik/(e-^^) 



bulk 



((1 



bulk 



(2.7) 



since the bulk average of a vanishes. As another example consider the concentric geometry, 
which may be viewed as a critical system inside a full spherical container (fc) with Hamilto- 
nian Tifc, where a hole S is generated by adding Tis- We are interested in the combination 
of free energies 



(-^5infc — -^fc) — (-^5 in bulk " -^1 



bulk J 



knTr, 



In (e~^s)fc + ln 



e-'^«)buik 



-keT, In (1 + or) 



f c : 



(2.8) 



the singular part of which is the function 5Fcon in ( |1.10|) . As a third example one may use 
( p.6a| ) also to study the behavior of (^E'(r)) in the concentric geometry as /2_ 0. One finds 



(^(r))c 



(^(r))fc {a ^(r))fc - (a)^ (^(r)) 



fc 



(2.9) 



to first order in a. This relation also applies for \1' Tnn, (see ( |2.3| )). Further, the expansion 
( p.6b|) holds also if the distant boundary is not a sphere concentric to S but, for example, a 
planar wall , and it should also apply slightly away from the critical point . 

Which quantities S appear in (|2.6b| ) depends upon the universality class of the surface 
of S. In case of the surface universality classes O and SB that do not break the 0{N) 
symmetry of the A^-component order parameter $ at the bulk critical point, all the H's are 
also 0{N) invariant, and the leading S is the energy density $^ For symmetry-breaking 
surfaces we only consider the case of a one-component order parameter $ and the class | or 
|. In this case S's, which are odd and even under $ —>■ — $ both appear in the expansion. 
The leading S is the order parameter $ and the leading even H is the energy density. 

In ( |2.6|) Rs is assumed to be large on a microscopic scale. The amplitudes A" for H = $ 
or $^ can be derived explicitly. They follow from the known result for the profile 
(\E'(r))5inbuik of the primary operators \E' = $ or $^ by comparing with (|2.7|) . From the 
leading behavior for |r — r^l ^ Rs, 



(Rs) 



(2.10) 



for \E' = $ or = respectively [jT4|. Here a denotes the surface universality class of S, 
A is the amplitude in the critical profile 



(^(ry,^)) 



(a) 

half space 



(2z)--* 



(2.11) 



of in a half-space with surface universality class a, and B is the amplitude of the bulk 
correlation function 



(^(r )^(0))buik = 5^r 



(2.12) 
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which also identifies the bulk exponent of 

Now we may calculate SF^on for small p as the singular part of (|2.8| ) by using ( p.6b| ) and 
( p.lOp , with the result 

6FaoM/kBTc -X^ ■ (^(0))fc = -^11^ p-* , (2.13a) 

where \1/ = $ if both surfaces a, b are symmetry-breaking and \1/ = $^ if one or both preserve 
the symmetry of the order parameter. Here we have identified Rs with R_ of Fig. la and 
have used the fact that the profile (\l/(0))fc at the center of a full spherical container with 
radius -R+ and surface universality class b is given by Af (i?+)~'^* (see [|19|)- Below we will 
also need the small-p limit 

Y{p) -2'^ S^^ X* p"* (2.13b) 

of the stress-tensor scaling function defined in ( |2.3| ). This follows either from ( p.l3a| ) by 

means of ( p.2|) or by inserting (|2.6b|) and ( p. 101) into with \E'(r) — T„„(ro) and using 
(TfcOfc = and 13 

(vl/(0) T„„(ro))fc = -ro"'^ (/?+)""* Sj' Af . (2.14) 

Eq. (P7[^ ) also implies, via (|rTIl| ) and (|r3|), the form - (2/t)-^* for the 

scaling function J-' in (|1.8| ) if k ^ 1. This leads to the Casimir interaction |P 



5* V -512 



5F(si2, Ri, R2) -ksT, (2.15a) 



between distant spheres (si2 ^ i?i,i?2) and 



<5Fspw(I^, i?i) ^ -fcsT, j (2.15b) 

between a planar wall and a single distant sphere {D ^ i?i). Here Eqs. ( p..4|) and (|1.9|) have 
been used. 



D. Density profiles and sliort-distance expansion near a planar surface 

In this section we study the profile (\l/(r')) as r' approaches one of the boundaries of 
the critical system. The normal component (T| 1 (r^)) of the stress-tensor average is 
regular near the boundaries ||38|. For the concentric geometry, this is apparent from ( p. 3D 



and it holds quite generally |3^,^. The behavior of the order parameter and the energy 



density is quite different, however, and is described by power-law exponents. In general 
these are determined by the scaling dimensions of surface operators and are different from 
bulk exponents Here we consider the particularly simple cases where \l/ approaches 

a ]9/anar boundary with universality class b and (^I/,&) is either ($,t)) ("^^jT)) ot ($^,0). 
Then the leading surface operator is T±± P,p|, pr3| , |30| , ^ , which leads to a surface exponent d 
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27| , |28[] . The corresponding short-distance expansion about the planar surface ||26| has the 
form |0|,|l],|9|,|l3l 



(2.16) 



The planar surface is located in r '-space at r[_ = 0. The factor in front of the square bracket 
is the critical half-space profile as in ( p.ll ) for surface universality class b. C^^ is a universal 
short- distance amplitude which generally depends on b [^. Its behavior is known in the 
Gaussian model with d > 2 and in the $^-theory near the upper critical dimension d = 4 
3l|,[l3|. The ellipses inside the square brackets in (|2.16| ) denote contributions of higher order 



in r^. The expansion ( |2.16| ) holds in thermal averages or correlation functions where besides 
the boundary at = and \E' other boundaries and operators may be involved. 

We are particularly interested in the small r_J_ behavior of (\l/(r||, r_[)spw in the sphere- 
near-planar-wall (SPW) geometry, which has been introduced as case (ii) just below Eq. ( |1.2|) 
and in Eq. ( |1.9| ). For the stress tensor 



d/2 



[D p + 2i?i)] 



r!^^ + D{D + 2Ri 



y{p) 



(2.17a) 



with 



P 



D + Ri- JD{D + 2Ri] 



/Ri 



(2.17b) 



and D as defined in ( p..9|) . Here we have used the general transformation formula [p!q ,p9|] 

(T^(r')) = HrT' (T„„(r))eoN (2.18a) 

with the local scale factor 



b{r') 



dv 



dr' 



-i/d 



R ■ r ' r 

1 + + 



f2 



(2.18b) 



i?2 4i?2 

Setting R = R^ leads to the SPW geometry on the left hand side of ( ^.18a| ) with 

D{D + 2i?i) = riin ^liv = i^R+ f , (2-19) 



which follows from ( [1.21 ). Using the form ( ^.3] ) for T„„(r) on the right hand side of ( ^.18a] ) 
and letting r _R+ so that — > leads to ( p. 17a]) . Eq. ( p.l7b|) follows e.g. from (p.. 51), 
since the quantity k from (|1.4| ) equals 1 + D/Ri in the SPW case. 

We close this section by introducing some of the tools needed in Sees. |ITl| and |V[ Critical 
systems with surfaces may be described [^| by a configuration probability ~ exp(— 7i) for 

the fluctuating A^-component order parameter $, where Ti separates into a bulk part 



bulk 



u 
4! 



(2.20) 
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and surface parts Tisurf- While the integration in (|2.20| ) is over the interior of the critical 
system, e.g. over R_ < r < i?+ in the concentric geometry, the surface contributions involve 
integrals (of powers of $ and derivatives thereof) over the surfaces of the system. Since we 
shall only consider systems right at the bulk critical point, we have not included a <l>^ term 
in H. 

We will also need the explicit form of the stress tensor at an arbitrary interior point of 
the system. It is given by 



(2.21a) 



with the canonical tensor supplemented by the improvement term 



4 d 



(2.21b) 



The quantity denoted by 0{u^) is needed to renormalize the stress tensor Since it is of 
third or higher order in the renormalized coupling ur (compare Eq. ( p.25| ) below), it does 
not appear in our explicit low-order calculations. 

The thermal average of the stress tensor in the concentric geometry takes the rotationally- 
invariant form 



d 



(2.22) 



Here r^^\ r^^^ are dimensionless quantities which only depend on the position via 



r . 



A fundamental property of the stress tensor is the continuity equation In the present 

case 







E 

I 



dr. 



-d-2 



dr 



(2) 



(2.23) 



which implies the important relation 



r — r'^(r„„(r))< 
dr 



(2.24) 



between the r-dependence of the radial component and the trace of the stress tensor. 
Eqs. ( |2.23| ) and ( p.24| ) hold for arbitrary m, rf, r and for arbitrary (rotationally-invariant) 
conditions at the surfaces. If the trace vanishes, as happens if T^buik is at the crit- 
ical fixed point of the renormalization group (RG), Eq. ( [2.24] ) implies that r'^(T„„)coN is 
independent of r. This is consistent with Eq. (|2.3| ). 

To set up the RG, we use reparametrization by minimal subtraction of poles in e 



45| , |46|p^ , |26[| . We note in particular the relation 



(2.25) 
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between the $^ coupling constant u and its renormalized counterpart ur. Here fi is the 
inverse length scale that determines the renormalization flow, and the (9('U/j)-contributions 
contain pole terms in e. The stress tensor ( p.21| ) is already renormalized and quantities 
such as )C in Eq. (|A8|) or the functions 



r 



(r/ R+,p; fir, ur) 



(2.26) 



with i = 1,2 from ( |2.22| ) are finite for e ^ order by order in ur. Eq. ( |2.26| ) follows from 
dimensional analysis, assuming that the boundary condition a, b take their fixed-point form 
and do not introduce any new lengths. If, in addition, the coupling ur equals its critical 
fixed-point value 



u 



R 



N + 2 
N + 8 



(2.27) 



the renormalization-group equation ||^ implies that the r^*) are independent of /ir. By Eqs. 
( p.22| )- (|2.24| ) and the vanishing trace, the r*^*^ are also independent of r/R^ and i. Thus, we 
find that 



g« {r/R^,p;f,r,u},) = g*{p) 



(2.28) 



only depends on p, d and the surface universality classes a, h. This provides another deriva- 
tion of (pi 



III. SITUATIONS WITH BROKEN SYMMETRY: MEAN FIELD 

APPROXIMATION 

In this section we consider the surface universality classes (a, h) = (T, T), (T,!), ihO) 
and {^,SB), where the order parameter symmetry is broken at the bulk critical point. With 
applications to the liquid-gas or consolute critical points in mind, we only consider a one- 
component order parameter. In the standard approach, the order parameter 

$(r) = (<l>(r)) + ^(r) (3.1) 

is decomposed into a non-vanishing average ($) and the fluctuation (p around it. Then a 
fluctuation (or loop) expansion is carried out, which amounts to an expansion in powers of 
u. 

To leading order one may neglect fluctuations altogether and determine ($) = ($)'''^^ by 
minimizing Ti. and ignore if. In the concentric geometry ($) only depends on the distance r 
from the center and [pl[] 



d — 1 n 

m H m = m , (3-2) 



with 



m = yw/e ($)(°) (3.3) 
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which has the dimension of an inverse length. For the trace of the stress tensor one finds 



k 

while the radial component is given by 



mTkk{v))^^l = -\,m'/u, (3.4) 



where 



r'^(TUr))(l = ^VM, (3.5) 



{2m^ - m^) + 2{d-2)mmr^ . (3.6) 



Here m = dm/dr, and the superscript zero denotes the 1/u contribution to (T), which is 
obtained by neglecting fluctuations. The first and second term on the right hand side of 
( ^.61 ) come from the canonical tensor and the improvement term, respectively. 

The quantities ( |3.4|) and ( p.5| ) determine the 1/u contributions to the right and left hand 
side of Eq. ( |2.24| ), respectively, and allow one to explicitly check ( |2.24| ) to this order. Note 
that for e — with -u > fixed the trace vanishes. Thus J4 is independent of r and represents 
a first integral of the differential equation ( p.2|) if d = 4. 

This integral was encountered in Ref. where it is pointed out that Eq. ( p.2|) is a 
special case of the generalized Emden-Fowler (GEF) equation For d = 4, (|3.2|) belongs 
to a subclass of the GEF equation, for which the first integrals are known explicitly ^8],^ 



A first integral for the mean-field equation in the spherically symmetric geometry may 
also be constructed by means of the stress tensor if in (2.20) is replaced by an arbitrary 
monomial interaction The leading contribution in u of the stress tensor always has 

a vanishing trace at the upper critical dimension d = 2M/{M — 2) (see Ref. |0). AH the 
corresponding mean-field equations belong to the above-mentioned subclass. 

Now we assume that the boundary conditions for m at r = i?_ and r = i?+ take their 
fixed-point forms. Then the dimensionless quantity Id in (|3.6| ) can only depend upon two 
ratios, say 

/, = /,(r/i?+, p) (3.7) 

of the three lengths r, In particular, for d = A 

h = I{p) , (3.8) 

i.e., depends only on p. Note that the term in r'^ (Tnn) con of leading order in uji, which 
follows from substituting ( ^.7| ) and ( |2.25| ) into (P75|), does indeed have the form ( ^.2(j| ). 



The scaling function Y{p) in (|2.3|) follows from (T„„)con for u = u*^^. Substituting ( p. 251) , 
( p. 27] ) and ( |3.8| ) into ( p.5| ) and keeping only the leading e-dependence, one finds 

y{p) = ; ^ I{P) + 0{e') . (3.9) 

To obtain solutions m in ci = 4 with the first integral 1^ = I, we follow the procedure in 
Ref. [^. With the substitution m = v/r Eq. ( |3.6|) becomes 
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dlnr = ±^2 dv /W{v) , 



(3.10a) 



with 

W{y) = VI + 2v^ + f 4 , (3.10b) 
and solutions can be obtained in terms of elliptic integrals. 

Case (i): For (a, b) = {], t) the order parameter profile m tends to +oo near the two surfaces 
and has one minimum in between. This implies J < and a solution with the parametric 



form 



V = ^ (3.11a) 

cos if 

In ^ = 7~^/' F U [(1 + lh)/2]'/') , (3.11b) 
where F is the elliptic integral 

F{ifi, k) = d(^ (l - sinV) "^^^ , (3.12) 

and 

7 = Vl -/. (3.13) 

Here r = Rm is the radius with minimum v, corresponding to = 0. At ip = ±7r/2 where 
the profile diverges, r takes the value r = R±. Eq. ( p.llb|) implies R_R^ = R?j^ and 

lnp = -/i(/) (3.14a) 

with 

h{I) = 2r"'K[[{l + lh)/2f^) . (3.14b) 

Here p is from ( p..6|) , and is the complete elliptic integral 

K{k) = F{'K/2,k) . (3.15) 

Eqs. ( |3.14| ), ( p.l3| ), and ( |3.9| ) determine the scaling function Y{p) for the stress tensor in 
( p.3|) . Note that / varies from to — oo as p varies from to 1. 

Case (ii): The solution for (a, b) = {SB, |) or (|, SB) is readily obtained from the solution 
in (i). Identifying Rm in ( |3.11b| ) with the radius of the SB boundary, i.e. considering 

i?^^ = R_ <r <R+ (3.16a) 

or 

R^<r < Rm = R+ (3.16b) 
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leads to the relation 



Inp 



(3.17) 



which only differs from (|3.14a|) by a factor 1/2. The reason for the identifications (p.l6|) 
can be understood by conformally transforming the mean-field (MF) solution (|3.11|) for 
(t, T) ill concentric geometry to spheres of equal size (SES) in r'-space, with the mid-plane 
= corresponding to the sphere with r = Rm (see the remark just below Eq. ( |1.2D ). The 



transformed MF solution is invariant under r'j_ —r'j_, and the normal derivative vanishes at 

r_J_ = 0. Since the MF equation A($)(°^ = (u/Q) (^{^)^^^^ is conformally invariant in = 4, 
this is the correct MF solution for a half-space r'j_> in d = 4 with a planar SB surface at 
= and a spherical hole with f boundary conditions. Transforming this solution back to 
the concentric geometry leads to ( p.l6|) and ( |3.17|) . Note that for the two geometries (p.l6|) 
the radial derivative dlnm(r)/dr at the SB surfaces does not vanish but equals — l/i?_ 
and — l/i?+ for ( p.l6a| ) and ( p.l6b|) , respectively Thus the MF order-parameter profile 
increases (decreases) on approaching the convex (concave) SB surface from the interior of 
the critical system, which is reasonable. The MF values ±1 for the ratio of the logarithmic 
normal derivative and the inverse radius of curvature in d = 4 are consistent with results in 
Ref. m. 



Case (iii): For (a, b) = (|, t), the profile tends to ±oo for r R± and has a zero in between. 
This implies / > 0. Though the /-dependence of m{r,I) is perfectly regular at J = 1, it is 
advantageous to choose different parametric representations for the intervals < / < 1 and 
/ > 1, respectively. They are given by 



1-7)^/2 tamp 



with 



In^ 

Rq 



1/2 



1 + 7. 



F [2/(1 + 1/7)]^ 



/2 



for < / < 1 (with 7 from ( |3.13| )) and by 

V = I^^^ tan (fi , 

with 



Rq 



1/2 



(3.18a) 



(3.18b) 



(3.19a) 



(3.19b) 



for / > 1. The profile diverges at = ±7r/2, corresponding to r = and has a zero 
at V9 = 0, corresponding to r = _Ro = (R+R-)^^"^ ■ Eqs. ( P.18b|) and (|3.19b|) determine the 
relation p = p{I). A form which holds for both < / < 1 and / > 1 is 



lnp=-/3(/) 



(3.20a) 



with 
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Here 



1 + 



(l + 1-'/') /2 



1/2 



(3.20b) 



(3.21) 



Case (iv): The solution for {a,b) = (0,1) follows from that in case (iii) by identifying Rq 



with the radius of the O boundary, where the order parameter vanishes |2^ . This implies 



lnp = --/3(/) 



(3.22) 



The four functions p = p{I) in Eqs. ( 13T7| ), (lOoP , and QO^ ) determine via (gj) 

the scaling function Y{p) of the stress tensor and via (|1.8| ), (|1.1CI|) , and ( p.2|) the free energy 
and the Casimir forces for the four pairs (t, |), {SB, |), (|, |), (O, j) of boundary conditions. 
Integrating according to (|2.2|) yields the parametric form 



9S, 



7r^ e 



'if mill 

dJ 



(3.23) 



for the scaling function, with k = k{p) from (|1.5|) and p = p{I) defined as above. On varying 
p or K, JF, (T„„)coN; and J remain negative in the cases (i) and (ii) and positive in the cases 
(iii) and (iv), implying attractive and repulsive interactions, respectively. The remaining 
integration in ( p.23| ) has been carried out numerically. Results for the scaling function for 
the four combinations of boundary conditions are shown in Fig. 2. Since the prefactor in 
( ^.23| ) behaves like 1/e as c? — 4, we have plotted the finite quantity = \imd^4^ (4 — d) T . 

It is instructive to check that our solutions exhibit the limiting behavior discussed in 
Subsecs. [11 Ij| and [11 (J| . Consider first the limit p ^ 1, where |J| — > cxd. In this case 



(/i,/3)^(2|/rVS23/2/-V^)-ir(l/v^) 



(3.24) 



behavior ( [2.4b[ ) with d = 4 and that 



Substituting the corresponding I{p) into (|3l9|) , one finds that Y{p) does indeed have the 

K{l/V2)f (3.25a) 



1 3 



27r2 



and 



A5S| = A||/16, A|| = -4A^|, Aot = A||/16. (3.25b) 

The expressions in ( [3.25| ) reproduce the known results pi]JTB|] for the Casimir amplitudes in 
the parallel-plate geometry for the symmetry-breaking cases and d near 4. 
Next consider p —* 0, where |/| — 0. In this case 

(/i,/3)^-ln(|/|/64) (3.26a) 

which implies 
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64- 



-P, -P^ P, P^ 



(3.26b) 



in cases (i)-(iv). In order to compare with ( p^.l3b| ) we need the expressions 



' 2 ' '2 



d 



(3.27) 



for AA/B in our four cases. Here \E' = [$, $, as noted just below Eq. ( ^.133] ). The 
results ( p.27|) follow from the behavior 



± 
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uz 



2 ' 



UZ 



2 ' 



(3.28) 



to leading order in u of ($), ($^) in the half-space with a =t or | (see Eq. (|2.11|) ), from the 
corresponding behavior 



,2-d 



(3.29) 



of ($^) in the half-space with an a = O or SB surface, and from the leading-order contri- 
butions 



yd ' , 2 (5, r^"^ 
to the bulk correlation functions (cf. Eq. (|2.12|) ). Here 



(3.30) 



(3.31) 



and u must be replaced by ( |2.25| ) with the fixed-point value ( |2.27| ) with = 1 for ur. Using 
the results ( p.26b|) and ( p. 27 ), one finds by means of Eq. ( |3.9|) and from 
for (i — > 4 that our four solutions are indeed consistent with (|2.13b| ). 

Now we turn to the density profiles and their behavior near a boundary. Consider first 
the r —>■ i?+ behavior of m in the concentric geometry for d = 4 with universality class b =t 
for the boundary at -R+. This follows directly from the differential equation (|3.10a|) , which 
implies 



\n^ = V2 dv/W{v). 

T Jv 



(3.32) 



Here we have used f (-R+) = oo and assumed that v is larger than a possible minimum in 
v{r). For r close to v tends to infinity, and one may expand the integrand in (|3.32|) in 
powers of 1/v. The integration is elementary, and one finds 



1 la 

40 • 



(3.33) 



The first factor on the right hand side is the profile inside a full spherical container with 
radius -R+, which follows from ( p.32|) with / = 0. The second term in square brackets 
respresents the leading correction for r i?+ due to the inner spherical boundary a with 
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radius The prefactor la,-] denotes / in the above cases (i), (ii), (iii), or (iv), with a =t, 
SB, I, or O (compare Eq. ( p.26b|) ). The result ( 3.33 ) also follows from the parametric 
representations ( p.llb| ), ( |3.18b| ), ( p.l9b| ) for m{r). 

The corresponding expansion near the planar wall with universality class 6 =| in the SPW 
geometry introduced below ( |1.2| ) follows from substituting (|3.33|) into the transformation 
formula 



(r') = [6(r')] 



(3.34) 



for mean- field solutions in c? = 4. Here 6(r') is from ( p^.l8b| ) with R = This leads to the 
result 



rnsaPWt(r')Mpw| (r|) 



(3.35) 



with 



2 = 4r|/2+/ {r!,Y + {2R. 



for the leading behavior of — r(r')] / R^ 



as r 



(3.36) 



0. Here we use the relation between 



(fc) 

the profiles mpw in half-space and m\ inside a full spherical container that follows from 



the conformal transformation ( p.34| ) |1l9| . Using Eqs. (|2.17a| ), (|2.17b|) and (|3.9| ) for (Tj_|_)c 
and the result 



(3.37) 



for the universal short- distance amplitude, we have checked that Eqs. (|3.35|) , ( p.36|) are fully 
consistent with the short-distance expansion ( |2.16| ) when applied to a (\E')spw = ('^')spw 
average for d — >• 4. 

In the other cases such as (\E', b) = ($^, O) with a =|, the short distance expansion 
( p.l6| ) about the planar wall with universality class b can be checked in a similar way ||51|| . 
The short- distance coefficient turns out ||31[] to be of order e° (see Eq. ( |4.24| ) below). 
As a consequence, the first term in square brackets in ( |2.16|) leads to a contribution which 
is smaller by a factor e than the contribution from the second term and is not present in the 
mean-field solution. 

In Appendix |B 1| we evaluate density profiles for i?_ for arbitrary fixed < r < i?+ . 
For p — > 1, the density profile in the parallel-plate geometry [|ll|] is obtained. This has 
been shown for (a, b) = (t, t) in Ref. and similar arguments apply for other a, b. 



IV. SITUATIONS WITH SYMMETRY PRESERVED: GAUSSIAN MODEL 



The boundary conditions (a, b) = (O, O), (O, SB), {SB, SB) do not destroy the 
bulk 0{N) symmetry of the A^-component order parameter at the bulk critical point. In 
these cases the starting point of the d = 4 — e expansion is the critical Gaussian model. 
This is defined by Eq. ( ^201) with 



u 



0. It is instructive |^2[ to consider the critical 
Gaussian model for arbitrary d > 2. For all d it is at its critical fixed point, i.e. dilatation 



17 



and conformally invariant on all length scales. This also holds in the presence of boundaries 
a and b if the boundary conditions take their fixed-point forms. 

In the concentric geometry the Gaussian order-parameter correlation functions or 'prop- 
agators' with a, 6 e {O, SB}, 



($i(ri)<l>j(r2))coN = 5ij G ab{r i,r2) = 5ij Gab{ri,r2;'n) 



(4.1) 



only depend on the positions ri,r2 through their distances ri,r2 from the center and the 
cosine of the enclosed angle 



T] = ri- r2/(rir2) . 
Here i and j denote components of the order parameter. One finds 

oo 

Ga,{rur2;v) = S, {r< r^y' Y.^^ iv) <i Xa{{r</ R-)') Xi>{{r>/ R+Y 

with 



1=0 



and 



min(ri,r2) , r> = max(ri,r2) 



^ = (d-2)/2. 



The are hyperspherical (Gegenbauer) polynomials and 



and 



xo{y)=y--, xsB{y) = y + - 
y y 



(4.1a) 



(4.2) 



(4.2a) 



(4.2b) 



(4.2c) 



(4.2d) 



The quantities W depend on the ratio p = R-/R+ of the radii of concentric boundaries and 
are given by 



W, 



(0 



w, 



Ofi 

(0 

0,SB 

(0 

SB,0 



p-1 _ p 



p-^ + P 



W. 



ii) 

SB, SB 



P-^ + P 



P-^ -P 



(4.2e) 



with 



P = p\ 



(4.2f) 
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The Gafe(ri, remain unchanged for ri ^ r2 and satisfy Gab = —S{ri — with A 
and S the Laplacian and the Dirac 5-function, respectively, in d dimensions. Note that xo 
and dxsB/dy vanish at the boundaries where y = 1. This imphes the fixed-point boundary 
conditions 

Go,b = , Ga,o = (4.3) 

and ED 



dr^ GsB,b = —75— GsB,b , f^r> Ga,SB = Ga,SB ■ (4.4) 

Here the first and second relations in ( [4 .31) and ( [4 .41) hold for r< = i?_ and r> = 
respectively. 

A conformal transformation 

G?r^ (r(, r^) = [b K) b {r^r'> G^b (ri, r^) (4.5) 



to the sphere near the planar wall (SPW) geometry shows that ( [4.4|) is equivalent to the usual 
Neumann boundary condition for a planar SB surface: Consider e.g. ( [1 . 1| ) with R = R^, 
which transforms the outer concentric boundary b = SB to the plane = (see Fig. 1). To 
prove that the derivative vanishes at r[j_ = for r2_|_ > 0, we show that the continuation in 
r{_|_ of G^^^^^ across the boundary is even in r{^. The points corresponding to r{ = (r^y , r{^) 

and its mirror image r[^^^ = (r^y , —r[j_) are ri and 

r = n Rl/rl , (4.6) 
respectively, as follows from inversion about the sphere r = Since 

X5B((rfV^+)') = XsB{in/R+)') (4.7) 

and 

r['h{ri^'^) =nb{r[) , (4.8) 

Eq. (I^D shows that G^^^^^ is indeed even in r{_|_. 

An inversion about the sphere r = y/R_R^ exchanges the boundaries a and h and implies 

G^(rur.,,)=(^)' gJ^,^,,) . (4.9) 

V rir2 / V ri r2 / 



One easily verifies that Go,sb, Gsb,o in (^3) consistent with ([4.91). 

Gab in (|4.2| ) should be compared with the propagator in the unbounded bulk 

oo 

Gbuik(ri, r2) = S, In - r2\-''^ = S, (r< r>)-'' ^ ^/''^(r^) (r</r>)^ . (4.10) 

1=0 

This equation follows from Gab in ( [4.2|) in the limits -R_ 0, -R+ ^ oo with ri, r2 fixed. 
Below we will need the difference 
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Gab{ri, r2) - Gbuiklri, = Sd 

1=0 



-1 



+ 



iL ^ ^(2) (^1^2)' (3) 
^■&+\ ' '^ab j^2\ ' "ab 



Rf 



{rir2y 



(4.11) 



which is regular for ri 



r2 at interior points. Here 
^^'^-(1,-1,-1,1) 
(-1,-1,1,1) 



a 



(2) 



,(3) 
'ab 



(-1,1,-1,1) 



for 



ab = {0 0,0 SB, SB O, SB SB) . 



(4.11a) 



(4.11b) 



In order to calculate (T„„(r))coN in the Gaussian model, we use (p.21a|) with -u = in the 
Hamiltonian ( |2.2(]| ). Since the average {Tki)hnik in the unbounded bulk vanishes, we evaluate 
the averages ($$)con and ($^)con that arise from Eq. (|2.21a|) by using Gab — Cbuik from 
( [4.11| ) instead of Gab- This yields 



CON 



00 



1=0 



-1 /rf - 3 + 



21 {d-2 + l) + -{d-2f 



(4.12) 



Using (|2.2| ), one eventually obtains 



d-2 



Einfi 



(1) 

ab 



P(« 



,2A 



1 

2A 



^rf - 3 + / 



^ |2/(rf-2 + /) + ^(rf-2)2| 



(4.13) 



where the relation between k and p is given in (|1.5|). The signs of both and (T„„)con 
are opposite to the sign of This holds for all values of k or p. Thus the interaction is 
always attractive if a = 6 and repulsive if a 7^ 6. We have evaluated numerically the sum in 
( [4.13| ). The function —To,o = —^sb.sb for rf = 4 is shown in Fig. 2. 

The profile of the energy density |^ is obtained from Gab — Cbuik- It takes the form 



00 



1=0 



-1 /d-3 + t 



2 + 



r 

r1 



2A 



a 



(3) 
ab 



r 



2A 



(4.14) 



Note the different sources of r-dependence in the expression ( |4.14| ) for the energy den- 
sity. The r-dependence is very similar in each of the contributions to (T„„) coming from 
Eqs. ( |2.21| ). However, (r/_R_)^^ and {R+fr)"^^ cancel in the sum, and (r„„) in ( [4.12| ) turns 
out proportional to r~'^, in agreement with ( p.3|) . The ((9$)^ terms with d parallel or per- 
pendicular to r that arise from the canonical tensor in (|2.21a|) contribute half of the first 
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term in the curly brackets on the right hand side of ([4.12|) , while the second term comes 
from the improvement J in (|2.21b ). The sum of the two terms in curly brackets equals 2A^. 

From ( ^.12| ) and ( |4.14| ), one may again check that the two combinations (a, 6) = (O, SS) 
and [SB^ O) are related by an inversion about the sphere r = yJR^ This operation 
preserves (T„„), while in ($^) the ratios r/i?_ and R+/r are exchanged. Note that (T„„) 
is the same for (O, O) and {SB, SB), too. However, the latter identity is limited to the 
Gaussian model, as we argue below. 

We now consider the limiting behavior of the stress tensor in ( [4.12|) . According to 



the discussion in Subsec. [II B| , (T„„) should approach the parallel-plate result when p 
1 — (L/R^) tends to 1. In this case 



P-^ ^ exp (2AL/i?^ 



(4.15) 



The exponential decay of the first factor in ( f4.12| ) is very slow, and the sum in ( |4.12| ) may 
be replaced by an integral 



E 



di 



-1 



2/^-Vr(d-2) . 



(4.16) 



Here P~^ is from ( [4.15| ) with A replaced by /, and the second and third factors in the sum 
of ( 4.121) have been replaced by their limiting forms for / 3> 1. The integral in ( 4.16 ) can be 
expressed in terms of the Riemann zeta function C{d) ||5^. Comparison with Eqs. ( p.3| ) and 
( p.4b|) shows that (T„„) does indeed have the expected limiting form, with the well-known 
Casimir amplitudes | ]IU| , P^ 



Ao,o = Ass, SB = -A^ (4vr) 




and 



A 



0,SB 



(1 



A 



(4.17a) 



(4.17b) 



for the A^-component Gaussian model bounded by plates of type O or SB. 

In the other limit p — > the leading behavior of (|4.12| ) comes from the / = term in the 
sum, which is given by 



E 



(4.18) 



This should be compared with the amplitude ratios 

AfAf/B^. = a« N/2 



(4.19) 



for the combinations ( [4.11b| ) of a and b. These follow from Eqs. ( p.29| ) and ( |3.3CI| ). Comparing 
with (2^) and ( p.31|) , one obtains the expected limiting behavior ( 2.13b ) with a;$2 = d — 2. 

We conclude the discussion of (T„„)con in the limits p — > 1 and p with two remarks: 
(i) Only the canonical part of T^i, i.e. the first term in curly brackets in ( [4.12| ), contributes 



in the limit p — > 1. This is to be expected, since the average of the improvement term, which 
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contains a total derivative parallel to the plates, vanishes due to translational invariance in 
the parallel-plate geometry. It is less obvious that the leading behavior for p — is only 
due to the improvement term (second term in curly brackets in ( [4.12| )), as we have noted in 
Eq. ( [4.18|) . (ii) For a normal critical system the <l>'^- interaction of ( p.20|) must be included. 
Then (T„„,) and thus Y and JF for (O, O) and {SB, SB) are no longer the same. From 
( p.l3b|) there should already be a deviation in first order in e = 4 — d for p — > 0, since 



{Afn / A, 



1 + 2e 



N 



(4.20) 



However, for p 



1, i.e. in the parallel-plate geometry there is no deviation in this order 
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Now we turn to the profile ( [4.14| ) of the energy density. This is the counterpart of the 
discussion at the end of Sec. |T|. For r fixed and i?_ — >■ only the second term in curly 
brackets in Eq. (4.14) survives, and the concentric profile tends to 



2-d 



CON 



a, 



(3) q (3) V- 



ab 



1=0 



r 



2A 



i?4 



r 



-i.d-2) 



iV-^($2(r))f,, (4.21) 



which is the profile inside a full spherical container (fc) with radius -R+ and universality class 
h (see the corresponding general expression in Ref. [|19|). Below we will need the difference 



N- 



($^(r))co.-($^(r)) 



fc 



2~d 



1=0 



-1 /d - 3 + r 



2 + ^2' 



r 



2A 



+ 



r 



2A 



(4.22) 



Note that the second term in curly brackets is different from that in ( [4.141 ). 

i?+ with 6 = 0, a^^fj = — 1. In this case the 



Here we study the behavior for r 



expression in curly brackets in ( [4.22|) tends to — 4A^ [(-R+ — r)/i?+] , and the sum over / can 
be expressed in terms of (T„„) in (f4.12|) . Using (|2.3| ), this leads to 



($2(r)) 



CON 



+ 



(4.23) 



for r — > i?+ and b = O. This is the counterpart of (|3.33|) . As in Eqs. (|3.34|) — (|3.37|) , one may 
transform to an SPW geometry and identify the leading coefficient in the short- distance 
expansion (SDE) about the planar wall with b = O. Now exponents —{d — 2) and d instead 
of —1 and 4 appear in Eqs. ( p.34| ) and (|3.35| ), respectively, and 



MO) 



/{NSd 



(4.24) 



This result, obtained here from the S^PWo geometry with a = O or SB, is consistent with 
the result for S^PWo geometry for c? — > 4, discussed in the paragraph after Eq. ( |3.37|) . It 
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is also in agreement with Refs. P,^, where the SDE of $^ about an O surface has been 
apphed to the profile in the parallel-plate geometry and in the half-space at T > and to 



3111 Ci*?^ for the Gaussian model was 



^$2 



correlation functions in the half-space at T^. In Ref. 
called B. 

Eq. ( |4.22| ) can also be used to evaluate the density profiles for arbitrary fixed < r < i?+ 
as R- —>■ 0, and one may check the predictions of the small-sphere expansion. This is done 
in Appendix |B 2 . 

Finally, we mention that for p ^ 1 the profiles ( ^4.14|) of the energy density reduce to the 



Gaussian results in the parallel-plate geometry P,pT 
analogous to those near Eqs. ( [4.15D and (|4.16| ). 



This may be shown with arguments 



V. SUMMARY AND CONCLUDING REMARKS 

We have analyzed the Casimir interaction between two spherical particles in a fluid at 
the critical point with boundary conditions 

(a, b) = (T, T), (T, i), (T,0), it SB), (O, O), (O, SB), {SB, SB) . 

The Casimir interaction is completely determined by a universal function JF = JF^ f, of a single 
variable k (see Eq. ( p..4|) ), which is a combination of the three lengths that characterize the 
configuration of two spheres. 

We have discussed the limiting behavior for spheres that nearly touch {k 1) and 
spheres that are far apart from each other {k —>■ oo) and calculated the complete functions 
J^ab{f^) for the case where the spatial dimension d is near 4 |]55| (see Fig. 2). From these 
results and corresponding information for (i = 2 and (a, h) = (|, t), (t, |), (t,0), (O, O) 
n|, one can make fairly reliable estimates of the scaling functions J-'ab in d = 3. 

Since k and are unchanged under conformal coordinate transformations, it is sufficient 
to consider a critical system bounded by two concentric spheres. We have derived ^(k) by 
calculating the average of the stress-energy tensor in this highly symmetric geometry. 

Besides the global free energy local thermodynamic averages, e.g. profiles of the order 
parameter and the energy density in the presence of two spherical particles, are also of 
interest. These can also be derived from the concentric geometry and we have pre- 
sented explicit results for d near 4 and for various surface universality classes a, b. We have 
also demonstrated that the stress-tensor average determines the modification of the density 
profile near the surface of one sphere due to the other distant sphere [P],pD|. 

The Casimir interaction and the density profiles for parallel plates are obtained as a 
special case of the results for concentric spheres when the radii -R_,-R+ tend to infinity 
with fixed difference i?+ — -R_ (compare Fig. la). The other extreme -R_ with i?+ 
fixed is also of interest. It determines the Casimir effect and the density profiles when one 
sphere is much smaller than other lengths in the problem. We have checked that the free 
energy contribution and the density profiles arising from the presence of the small sphere 
are consistent with a 'small-sphere operator expansion'. This method has been introduced 



in Ref. and is also explained in Sec. |II C| of this paper. The density profiles show still 



non-trivial scaling forms in this limit (cf. Appendix B). 
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We have also found simple expressions [^] for the interaction between distant spheres in 



an infinite host fluid if the latter is not at the critical point. Here both the host correlation 
length ^ and the distance s between the spheres have to be much larger than the radii of 
the spheres, but the ratio ^/s is unrestricted. Similar expressions have been found for the 
interaction of a single sphere with a distant planar boundary. 

We conclude by discussing possible experimental observations of Casimir interactions 
between spheres at the critical point. The case of fluid mixtures at their consolute point is 
of particular interest, since the attractive interaction between two immersed equal spheres 
decays extremely slowly, with a power law 

-kB^Aisu/Ri)-'-''^ (5.1) 

in (i = 3 as the separation Si2 becomes much larger than the radius Ri of the spheres. The 
exponent follows from Eq. ( ^.ISal ), with 2a;$ = 2/5/z/ = 1.04 in d = 3 [0]. In the only 



previous theoretical work we know on this topic, de Gennes obtained an extremely good 
estimate 1 for this exponent by minimizing a phenomenological free energy functional. The 

universal amplitude A = (^f) /-B* is estimated to be somewhat larger than \/2 [0]. The 
Casimir interaction ( |5.1| ) should be compared with the van der Waals interaction which 
decays much faster with an exponent —6 (non-retarded) or —7 (retarded). The prefactors 
are comparable. For example, for fused silica spheres in water the van der Waals prefactor 
is about —ks ■ 10^ Kelvin We note that for spheres in "^He at the A-point {XY model) 
the Casimir exponent is 2a;$2 = 2{d — l/u) = 3.02 for d = 3. Thus the decay is faster than 



for fluid mixtures (see (|5.1|)). 

Finally we mention two classes of experiments one could conceivably carry out to study 
the Casimir interaction: 

(i) Interaction between only two objects, for example a single sphere near a planar wall: By 
measuring the force on the sphere for varying distance from the wall, one can in principle 
obtain the whole scaling function T{k). Of course, the distance must not exceed the correla- 
tion length of the nearly critical host binary fluid. Correlation lengths of several 1000 A are 
feasible, and one would like to have a sphere diameter of a few hundred A. Measurements 
of the interaction potential of isolated pairs of spheres of somewhat larger size have been 



reported for a charge stabilized colloid pB|, see also 



(ii) Aggregation phenomena in colloids: The Casimir interaction between more than two 
spherical particles at large [^] separations is a sum of pair interactions, which have the 
gravitation-like form ( |5.1| ) for separations s smaller than the correlation length ^ of the pure 
host and vanish for s ^ ^. On approaching the critical point ^ strongly exceeds the radii and 
mean separation of colloid particles. This offers a possible mechanism for their aggregation 
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Note added in proof : After this work was completed we realized that the widely 
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used non-retarded Hamaker approximation 5E for the van der Waals interaction of two 
spheres is also conformally invariant (hke the critical Casimir energy 5F studied in this pa- 
per). This follows immediately from its representation — (A/vr^) / d^ri d^r2 |ri — Y2\~^ 
with integrations over the two spheres and allows to express 5E as a function of the 
invariant quantity k defined in (1.4). The explicit result [57] 6E {su, Ri, R2) = 
— {A/6) {{k, — 1)^^ + (k + l)"""^ + In [(k — !)/(/€ + 1)]} may be used for a first comparison 
of critical Casimir and van der Waals interaction. Note that for close spheres {k near 1) 6E 
has the same dependence on R2, D as 5F and can be obtained from the right hand side 
of (2.5) by replacing fc^Tc 5*^ ^ab by —A/?) and setting d = 



APPENDIX A: FREE ENERGY AND STRESS-TENSOR AVERAGE 



An infinitesimal coordinate transformation 

X = X + a(x) 

changes a geometry G to a geometry G. We consider two examples: 

(i) For the geometry of non-overlapping spheres in Fig. 1 b where x = r', 

ai{r) = a 6i^± 6 {r[_) , 



(Al) 



(A2) 



with the unit step function, implies a rigid shift of the right sphere in the _L direction and 
increases the distance S12 between the spheres by the amount a. 

(ii) For the concentric geometry in Fig. 1 a where x = r, we consider the transformation 

nn 



The corresponding change in the singular part of the free energy is given by ||T8 
Fq = Fg- ksT, J d'x J2 {dak{^)/dxi) {Tm{^))g 



(A3) 



to first order in a. Here Tki is the stress tensor. In case (i) 
d 



ds 



12 



SF{s,2, Ri,R2)/kBT, = J d"-V|| (r^(r|(,0)) 



with the average taken in the geometry of Fig. 1 b. In case (ii 
dak/dri = a 



dki&ir - To) + — 6{r - tq) 
r 



(A4) 



(A5) 



The result simplifies at the critical fixed point since the trace of the stress tensor vanishes 
for R_ < r < R+ JlSl; a possible boundary contribution at r = _R+ drops out in the change 
of 5FcoN [2^- Then only the second term on the right hand side of (KE) contributes and 
picks out the radial component T„„ of the stress tensor. This leads to 



— FcoN(i?-/[(l + «)i?4 
da 



-kBTc SdTQ {Tnn{ro)) 



(A6) 



a=0 
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to first order in a, which implies Eq. ( |2.2| ). 

Since both 6F{si2, Ri, R2) and 5Fcon{p) can be expressed in terms of the same function 
(see Sec. |), the ratio between the left-hand sides of Eqs. ( |A4| ) and ( pl^ ) is simply 
— ((9K/5si2)/(pdK(p)/dp), with K and from (|1.4|) and (|1.5|) , respectively. It is instructive 
to verify from the transformation formula p.l8a| ) for the stress-tensor average at the critical 
fixed point that the right-hand sides of ([A4|) and ( p.2|) lead to the same simple ratio. With 

= R ■ r' = and r = i? in Eq. ( |2.18a| ), an easy ry-integration leads to the ratio l/(2i?). 
This equals the above-mentioned ratio, since 



2R 



S12 
Ri R2 



p 



(AT) 



as follows from writing S12, -Ri, R2 in terms of the intersection points I, II, III, IV in (|1.2| ). 

Eq. ( |A^ ) is not restricted to the critical fixed point but applies for arbitrary values ur 
of the renormalized $^ coupling strength (see Eq. ( |2.25| )). This can be used to show the 



universality (Mij-independence) of the singular part of the free energy in the long-distance 
limit. For simplicity consider a situation where the boundary conditions have their fixed- 



point form. Then by naive dimensional analysis the right hand side of 

as 



(?1 



0)) = — /C(i?i/Si2 

Sl2 



R2/S12 ; /iSi2, Ur) 



can be written 



(A8) 



In the 'shift operator' / d'^^^rij 
contribute, and the average (^ 
malization group equation 

Pw- + (3{ur) 

O jJ, 



T\ I (r|j, 0) the improvement term J in (|2.21b| ) does not 
]) is renormalized and obeys a simple homogeneous renor- 



d 



duR 



/C (-R1/S12, -R2/S12 ; /iSi2, Ur) = 



(A9) 



with the usual /3-function [^,@ . Standard arguments then imply that in the long-distance 
(infrared) limit /iSi2 — > C)0 with Ri/si2 and R2/S12 fixed 



/C (-R1/S12, R2/ S12 ] pSi2, Ur) > /C (-R1/S12, R2/ S12; 1, 



u 



R) 



(AlO) 



with from ( |2.27| ). This shows the -u/j- independence of the asymptotic critical Casimir 
force (ptil). 



APPENDIX B: DENSITY PROFILES AND SMALL- SPHERE EXPANSION 



Here we study the density profiles (\E')con = {^)a,b in the concentric geometry in the 
limit R- with f = r/i?_|_ fixed at an arbitrary positive value < f < 1. We show that 
this form is indeed in agreement with |Q, i.e. 



(^(r)),,,-(^(r)) 



fc 



ijq/ 



(^(0)^(r))fe-(^(0))fc(^(r))fc 



(Bl) 



which follows from the small-sphere expansion (SSE) (p.9|), ( |2.6b| ), and ( p.lOp . As before 
the full container (fc) has a surface with universality class b. In ( pi|) , \1/ = $ if both and 
the correlation in square brackets are non-vanishing, and \1/ = <l>^ otherwise. 
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1. Cases with broken symmetry for d — > 4 



We start with the two cases (a 6) 
breaks the symmetry. Consider first ^ 



(1,0) and {^,SB), where only the inner sphere 
$. Then \l/ = $, and (|B|) imphes 



(0) _ /6 v{r) 



(B2a) 



with 



vir) 



2V2 



where 



=p r 



(B2b) 



Here the upper sign holds for 6 = O and the lower sign for h = SB. To obtain (p2|) we have 
used the vanishing of for b = O and SB and the result for ($(0)$)^ that follows via 
a simple conformal transformation |jl9[ from the half-space correlation function, which for 



c? — i> 4 is that of the Gaussian model with Dirichlet or Neumann boundary conditions p6 



We have also used ( p.28| ) and ( p.30| ) to evaluate Af and 5$, respectively. 



The predictions (|B^ ) of the SSE are consistent with the solutions of the differential 
equation (|3.10a| ). For b = O this has the form 



r Jo 



dv 
W{v) 



(B3) 



Note that / is the same for (a, b) = (t, O) and (O, |) and equals 64p^ by (|3.26b|) for our case 
p 0. It follows from ( |B2b| ) that the t;^-term in the function W in ( |B3D is of order and 



can be neglected, since I and 2v'^ are both of order p^. The integral in ( P3|) is then trivial 
and leads to the result ( [B2b| ) with the upper sign. 

For b = SB, the solution follows from ( [B3| ) on replacing the lower limit of integration by 
the positive w-value 



1/2 



1 



1/2 



(B4) 



for which W{v) vanishes. This guarantees that dv/dr vanishes at r = |61|]. For p — > 0, 
in which case |/| —>■ 64p^ by Eq. (|3.26b| ), one may again neglect inside W{v), and the 
lower boundary of integration becomes (|/|/2)^/^. The integral is simple and leads to the 
result (|B2b| ) with the lower sign. 

The leading contribution to the energy v]/ = $2 _j. 4 jg given by [B^ 



(B5) 



and follows immediately from (|B2D . Since ^^(O) = $^(0), the right hand side of is, in 
leading order, the square of its value in the previous cases, and Eq. ( [BlD is again fulfilled. 

Now we turn to the cases where the outer sphere breaks the symmetry, e.g. b =]". For 
the order-parameter profiles \I' = $ Eq. ([BI|) implies 



1 6 v\''\r) + Sva,^{r) 
u r 



(B6a) 
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with 



2^/2 



(B6b) 



and 



16^2 



q- + 



12 



1 _ ^ in(f-i) 



(B6c) 



with / from ( ^.26b| ) for the four cases a = (t, SB, |, O). Here ( |B6b| ) represents the non- 
vanishing profile ($)fc in a full container (fc). It can be obtained by means of a conformal 
transformation from the half-space profile ( |3.28| ) as described e.g. in Ref. \T9\. Eq. ( |B6cD 
respresents the right hand side of ([B 1|) . Here ^'(0) is ($, $, in the four cases (see 
remark just below ( p.27|) ). For d —>■ 4 



(<l>(0)<l>(r))fc-(<|.(0))fc($(r)) 



fc 



(B7a) 



is determined by Gaussian fluctuations (p about the mean-field profile The quantity 

Kf^ can be obtained from its counterpart in the half-space whose form is known (see Eq. 
(4.96) in Ref. The corresponding correlation for \E'(0) = $^(0) follows from 



($2(0) <|.(r))fc - ($2(0) )fc($(r))fc ^ 2 ($(0))f (r) 



(B7b) 



Note the limits 



and 



SVa- 



16^2 r 



7-. 



(B8a) 



6v, 



(fc) 



40 



'R4 



R4 



(BSb) 



which follow from Eqs. (|B^ ). Eq. (|B8b|) is consistent with the expansion (|3.33|) about the 
surface, which applies for arbitrary R_/R^. 

Again the prediction ( P6|) of the SSE may be checked by comparing with the solution 
of the differential equation ( |3.10a| ), which is now given by Eq. ( p.32| ). Since v tends to the 
finite positive value wj^^'' in ( P6b| ) as / ^ 0, one may expand the integrand in ( p.32|) with 
respect to / and invert the function r{v) to obtain v = v{r), order by order in /. In first 
order the result reads 



vir] 



V2 g- 



1-^ — 



(r) 



[Wo{v)\ 



(B9) 



where Wq is W for / = 0. The integral on the right hand side is straightforward, and one 
sees that Eq. ( p9D indeed implies the result (|B6|) . 
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2. Cases with symmetry preserved in the Gaussian model 



Here we consider the i?„ —>■ behavior of the left hand side of ( [BID ^ — '^'^(r) and 
a, be {O, SB}. The result 



^-2(d-2) 



1 + 



(BIO) 



follows for the Gaussian model from the / = term in ([4.23| ), 



that cr^^ only depends on a, while o"^'^'' only depends on h. This is in agreement with the 
right hand side of ([Bll). It is easily verified that ( |B10[ ) equals the right hand side of ( |Bip : 
The leading operator is \E' = and the $^ ■ $^ correlation function inside a full spherical 
container with an O or SB surface, which appears on the right hand side of ( |B1| ) in square 
brackets, follows from its counterpart in the half-space with Dirichlet or Neumann boundary 



.(3) 



(1) 

using al^; 



■ a 



(3) 
ab 



a 



(2) 
ab ■ 



Note 



conditions, respectively. The prefactor A* /B^2 in (pl|) is obtained from (|3.29|) and ( p.30D 
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form ( p.6|) , the leading S being a line integral dt\&(et) with the unit vector e in the 
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governed by an exponent larger than the bulk exponent d — 2, and the expansion does 
not apply.) In this case = and the amplitude A'/i?^*^-! = C{d) can be explicitely 
determined. For the interaction free energy 6F12 of two crossed wires with O (ordinary) 
boundary conditions, diameters 2Ri, 2R2, and a minimal separation su ^ Ri , R2, one 
finds 5Fi2 = —ksTc {R1R2/ s\2)^ C^V from the arguments which led to Eq. (15) in Ref. 
TA\. Here E = x<j,2 — l = c/— 1 — l/z/ = 0.51 for d = 3, N = 2, which should apply to 



^He. The quantity V = J^^ dh dts ($^(eti) $^(si2 + eta)) bulk ■ s?f is independent of Su 
but depends on the angle between the wires. 

[35] For simplicity consider T>Tc. Both the finite bulk correlation length ^ ~ (T — T^)^*^ and 
the 'large' distances mentioned right after Eq. ( 2.6b| ) have to be much larger than Rs in 



order that Eqs. ( p.6| ) apply. The amplitudes X^^, are independent of T — and are given 
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Figure Captions 



Fig. 1: Conformal mapping of concentric spheres with center at the origin of r-space (Fig. 
1 a) onto spheres in r '-space with centers on the hne r' = R/i? (Fig. 1 b) as imphed by 
Eq. ( IlTI) . The sphere denoted by the broken hne with radius R is mapped onto the plane 
r[_ = 0. The two points r = — R and r = R of this sphere are mapped onto r ' = and 
r ' = oo, respectively. The two spheres with radii -R_ and -R+ are mapped on spheres with 
radii Ri and R2, respectively. The corresponding intersection points of the two spheres with 
the coordinate axes parallel R are denoted by I, II, III, IV. For the case R- < R < R+ 
shown in Fig. 1 a the images r' = 2 R and r' = —2 R of r = and r = 00 are located inside 
the two spheres in Fig. 1 b and are closer to the origin than the centers of the spheres. 

Fig. 2: Scaling functions JF = J^^b in a critical Ising-type (A^ = 1) system for spatial 
dimension d —>■ 4. The function J-'osb, which is not shown, would nearly coincide with 
—^00 = —^SBSB in this plot. For the four symmetry-breaking combinations {ah) = (| |), 
(tt), (t O), (t SB), the functions behave as 1/e for d ^ 4, and we have plotted the 
finite limit ^ab = limcz^4(4 — d) Tab- 
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Fig. 1 
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Fig. 2 
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